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Here we examine classical 0{n) spin systems with arbitrary two spin interactions (of unspecified 
range) within a general framework. We shall focus on translationally invariant interactions. In the 
this case, we determine the ground states of the 0(n > 2) systems. We further illustrate how one 
may establish Peierls bounds for many Ising systems with long range interactions. We study the 
effect of thermal fluctuations on the ground states and derive the corresponding fluctuation integrals. 
The study of the thermal fluctuation spectra will lead us to discover a very interesting odd-even n 
(coupling-decoupling) effect. We will prove a generalized Mermin-Wagner-Coleman theorem for all 
two dimensional systems (of arbitrary range) with analytic kernels in k space. We will show that 
many three dimensional systems have smectic like thermodynamics. We will examine the topology 
of the ground state manifolds for both translationally invariant and spin glass systems. We conclude 
with a discussion of 0{n) spin dynamics in the general case. 



I. INTRODUCTION 

In this article we aim to unveil some of general proper- 
ties of 0{n) spins systems having two-spin interactions. 
We shall mostly concern ourselves with translationally 
invariant systems. 



two dimensions with an analytic interaction kernel never 
develop spontaneous magnetization. At the end of the 



The outline is as follows: In section (III) we introduce 
frustrated toy models rich enough to illustrate certain of 
the general features that we aim to highlight. These mod- 
els will be employed for illustrative purposes only. The 
bored reader is encouraged to skip the somewhat verbose 
exposition and merely study the two Hamiltonians. 

In section([V) we discuss the ground states of Ising 
spin models and show what patterns one should expect 
in general. Once the ground states will be touched on, we 
will head on to show how Peierls bounds may be estab- 
lished for many systems having infinite range interactions 
if the ground states are simple. In section([v|) we shortly 
review mean field solutions of the general two spin Ising 
models. 



In section(VII) we prove that, sans special commen- 
surability effects, the ground states of all 0{n > 2) will 
typically have a spiral like structure. 

In section( VIII) we will make an exceedingly simple 
spin wave stiffness analysis to gauge the effect of thermal 
fluctuations on the various 0(n > 2) ground states. 

In section([X) we will discus thermal fluctuations 
within the framework of "soft-spin" XY model. We will 
see that the normalization constraint gives a Dirac like 
equation. In the aftermath, the fluctuation spectrum will 



be seen to match with that derived in section( VIII). We 
will show possible links to smectic like behavior in three 
dimensions. 

Next, we go one step further to study the fully con- 
strained "hard-spin" 0{2) and 0(3) models and show (in 
section ^ that all translationally invariant systems in 



section our analysis will match that of sections (VIII) and 

H) 

We extend the Mermin-Wagner-Coleman theorem to 
all two dimensional interactions with an analytic inter- 
action kernel in momentum space 

In section(XI) we will examine the "soft-spin" version 
of Heisenberg spins. We will see that it might be naively 
expected that the spin fluctuations in odd n spin systems 
are larger than in those with an even number of spin 
components. One of the fluctuating spin components will 
remain unpaired. 



Next, in section(XII), we carry out the spin fluctuation 
analysis for four component soft spins to see that their 
spectra coincides with that predicted in the earlier spin 
stiffness analysis. 

In section(XIV) we compute the critical temperature 
of all translationally invariant 0{n > 2) spin models with 
mean field theory. 



In section(XIII) we show that in the limit of large n 
both odd and even component spin systems behave in 
the same manner. Essentially, they all tend towards an 
"odd" behavior. 



In section(XV) we briefly remark that much of anal- 
ysis is not changed for arbitrary two spin interactions 
(including spin glass models). 

We conclude with a discussion of 0{n) spin dynamics. 

A central theme which will be repeatedly touched on 
throughout the paper is the possibilities of non-trivial 
ground state manifolds. If the system is degenerate the 
effective topology of the low temperature phase of the 
system may be classified in momentum (or other basis). 
In such instances the low temperature behavior of the 
systems will be exceedingly rich. 
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II. DEFINITIONS 

We will consider simple classical spin models of the 
type 



(1) 



Here, the sites x and y lie on a (generally hypercubic) 
lattice of size N. The spins {S{x)} are normalized and 
have n components: 



(2) 



at all lattice sites x. 

We shall, for the most part, consider translationally 
invariant interactions V{x,y) = V{x — y). 

We employ the non-symmetrical Fourier basis conven- 
tion 



(/(fc) - E.-^(a 



~ik-x . 



F{x) 



wherein the Hamiltonian is diagonal and reads 



where v{k) and S{k) arc the Fourier transforms of V{x) 
and S{x). 

More generally, for some of the properties that we will 
illustrate, one could consider any arbitrary real two spin 
interactions (x|V^|j7) which would be diagonalized in an- 
other basis {\u)} instead of the Fourier basis. 

For simplicity, we will set the lattice constant to unity- 
i.e. on a hypercubic lattice (of side L) with periodic 
boundary conditions the wave-vector components fc; — 
where r/ is an integer (and the real space coordinates 
xi are integers). 

Throughout this work we will use K to denote recipro- 
cal lattice vectors and A(/c) as a shorthand for the lattice 
lattice Laplacian: 



A(fc) = - cosfci 



(4) 



1=1 



In some of the frustrated systems that we will soon 
consider, v{k) may be written explicitly as the sum of 
several terms: those favoring homogeneous states k 
0, and those favoring zero wavelength k —^ oo (or k —i- 
(tt, TT, tt) on a lattice.) As a result of this competition, 
modulated structures arise on an intermediate scale. 



III. TOY MODELS- FOR ILLUSTRATIVE 
PURPOSES ONLY 

Although we will keep the discussion very general, it 
might be useful to have a few explicit applications in 
mind. There is a lot of physical intuition which underlies 
the upcoming models. Unfortunately, insofar as we are 
concerned, they will merely serve as nontrivial toy mod- 
els on which we will able to exercise our newly gained 
intuition. 

The systems to be presented are frustrated: not all two 
spin interactions can be simultaneously satisfied. 

We choose these rather nontrivial toy examples as they 
highlight some possible richness which is typically absent 
in the more standard spin models. As such, they will 
point to typically ignored subtleties. 



The Coulomb Frustrated Ferromagnet 

Let us introduce our first toy model "The Coulomb 
Frustrated Ferromagnet" . 

This is a toy model of a doped Mott insulator, where 
the tendency, of holes, to phase separate at low doping 
is frustrated, in part, by electrostatic repulsion In 
three dimensions, a simple spin Hamiltonian which 
represents these competing interactions is 



Mott 



S{x)S{y) 



Stt la; - y\ 



= ^E[^(^)+Ei^-^i"']|^(^)i' 

r. K 



(5) 



Here, S{x) is a coarse grained scalar variable which rep- 
resents the local density of mobile holes. Each site x 
represents a small region of space in which S{x) > 0, and 
S{x) < correspond to hole-rich and hole-poor phases 
respectively. In this Hamiltonian, the first "ferromag- 
netic" term represents the short-range (nearest-neighbor) 
tendency of the holes to phase-separate and form a hole- 
rich "metallic" phase, whereas the frustrating effect of 
the electrostatic repulsion between holes is present in the 
second term. Non-linear terms in the full Hamiltonian 
typically fix the locally preferred values of S{x). One 
may consider d ^ S dimensional variants wherein the 
spins lie on a hypercubic lattice, and the Coulomb kernel 



in Ho is replaced by - y\'^ (or by In \x-y\] 

in two dimensions) where Vld — /T{d/'l). Here the 
competition between both terms, when Q ^ 1 favors 
states with wave-numbers ~ Q^^^. The introduction of 
the Coulomb interaction is brute force non-perturbative: 
it is long range. Moreover, the previous ferromagnetic 
ground state becomes, tout a' coup, infinite in energy. 
We will, for the large part, focus on the continuum limit 
of this Hamiltonian where the kernel becomes 
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VMTuik) = Qk-' + k'[l -^J] (6) 

After rescaling, this may also be regarded as the smaU k 
limit of the more general 

VQ{k) = A(fc) + Q[A(fc)]-i + A[A{k)]^ 

+A^(1 -cosfc,)(l -cos/cj) + 0(fc6). (7) 

The constants A and A are pinned down if we identify 
'^q{^) = VMott{k)- Here, we will modify them in order 
to streamline the quintessential physics of this system. 
First, we set A = Q: In the continuum limit this term 
is not large nor does it lift the "cubic rotational symme- 
try" of the lattice (i.e. those transformations which leave 
A(fc) invariant) present to lower order. Next, we allow 
A to vary in order to turn on and off "cubic rotational 
symmetry" breaking effects. 

Note that VQ{k) with may be regarded as vmou aug- 
mented by all possible next to nearest neighbor interac- 
tions. As our Hamiltonian respects the hypercubic d ~ 2 
point symmetry group, by surveying all possible values of 
A we should be able to make general statements regarding 
the possible phases (within the planes) of real doped Mott 
insulators. When A > the minimizing wave-vectors will 
lie along the cubic axis and "horizontal" order will be 
expected. When A < the minimizing wave-vectors lie 
along the principal diagonals and diagonal order is ex- 
pected. At large values of Q, when the continuum limit 
no longer applies, trivial extensions of these minimizing 
modes are encountered where one or more of the wave 
vector components is set to tt. 

Note that if the ferromagnetic system were frustrated 
by a general long range kernel of the form V{\x — y\) ~ 
\x — y\~P we could replace the [A(^)]~^] in VQ{k) by the 
more general [A(fc)](^'~'')/^. Here, in the continuum limit, 
the minimizing modes are ~ Q'^/i^+'^-p) and as the reader 
will later be able easily verify all our upcoming analysis 
can be reproduced for any generic long range frustrating 
interactions with identical conclusions. 

Shown above is the manifold (M) of the minimizing 
modes in k space. When no symmetry breaking terms 
(A — 0) are present, in the continuum limit it is M is the 
surface of sphere of radius Q^^^. If A ^ this degeneracy 
will be lifted: only a finite number of modes will minimize 
the energy. When A > there will be 2d minimizing 
modes (denoted by the big X in the figure) along the 
coordinate axes. In the up and coming we will focus 
mainly on A > 0. When A < 0, a moment's reflection 
reveals that there will be 2** minimizing modes along the 
diagonals, i.e. parallel to (±1,±1,±1) (and in this case, 
they will have a modulus which differs from Q^/**). 

Unless explicitly stated otherwise, we will set A = for 
calculational convenience and when a finite A is invoked 




FIG. L The shell of modes which minimize the energy in 
the continuum limit. For the case just discussed This sphere 
is of radius |g| = Q^^"^. 



it will be made positive (to avoid the A dependence of \q\ 
incurred when the former is negative). At times, we will 
present results for VQ{k) at sizable Q, even though the 
model was motivated as a good caricature of vmou only 
in the continuum limit (at small wave-vectors ~ Q^^^)- 
Membranes 

In several fluctuating membrane systems, the afhnity 
of the molecular constituents (say A and B) for regions of 
different local curvature frustrates phase separation 

Let us define S{x) to be the difference between the A 
and B densities at x. 

In the continuum, the energy of the system contains a 
contribution, 

H^^^ d^x |V5|2 (8) 

reflecting the demixing of A and B species. Instead of 
considering long-range interactions, we now allow for out- 
of-plane (bending) distortions of the sheet. Specifically, 
we assume that the two molecular constituents display 
an affinity for regions of different local curvature of the 
sheet. This tendency can be modeled by introducing a 
coupling term between the local composition S{x) and 
the curvature of the sheet. 

Provided that the distortions remain small, we may 
write 

Hc= [ d?x [i(7|V/i(f)|2 + ^[v2/i(x)]2 + A5(f)V2/i(f)] 



3 



where h{x) represents the height profile of the sheet (rel- 
ative to a flat reference state), a is its surface tension , 
and K, is its bending modulus; A, the coefficient of the 
last term in the expression measures the strength of the 
coupling of the local curvature V^h and the local compo- 
sition (/), which wc have included here to lowest (bilinear) 
order. This coupling term reflects the different affinities 
of the molecular constituents A {S = 1 corresponds to 
pure A composition) and B (S = corresponds to pure 
B composition) for, respectively, convex (V^/i > 0) and 
concave (V^/i < 0) regions of the interface. We now min- 
imize the total energy H = H^ + Hc, w.r.t the membrane 
shape {h{x)}. 

the variational cqns, where in obtaining to the last line 
wc have employed S{dfh{x)) = dfSh{x), 5{dih{x)) = 
diSh{x), and integrated by parts twice. 
Thus 

-aV^h + KV^{V^h) + AV'^S = 0. (9) 

If |kV^(V^/i)| < min{ I o-V^/i I, \AV^S\} then an approxi- 
mate solution to the last eqn is 

AS w ah, (10) 

and in Tic, after an integration by parts, 

ASV^h « AS-V^S - — {VSf. 
a a 

Hmi. + Hc^J d'x [^6'|V5|2 + ^{V'Sf] 

A^ 

where b' = b . (12) 

(7 

This effective energy reads 

H = j (fk Vmembrane{k)\(p{k)\'^, (13) 

where Vmembraneik) = ^k"^ + f^fc** is the 2D Fourier 
transform. A negative b' obtained when b < A^/a, sig- 
nals the onset of a curvature instability of the sheet. This 
instability generates a pattern of domains that differ in 
composition as well as in local curvature and thus as- 
sume convex or concave shapes. The characteristic do- 
main size corresponds to the existence of the minimum 
of the free energy at a non-zero wave number. The mod- 
ulation length d ~ ^/{A^K/a'^)/\b'\. 



After scaling, this model may be regarded as the con- 
tinuum version of the frustrated short range kernel 

v^{k) = zA^{k)-A{k) (14) 

(where z = — A^ n / {a^h')) on the lattice. 
The real lattice Laplacian 

,^1 . I {2d for X — y ^, 
^^1^1^^- 1-1 for||f-,1loo = l ^''^ 

Notice 

that {x\A^\y) = for \\x - y\\oo > R {= Range) . Our 
system is of Range = 2. 
Explicitly 

{x\A^\y} = 2d{2d + 2) for x = y 
—Ad for |x — j7| = 1 

2 for {x — y) = {±ee ± eg') where £ ^ £' 
1 for a ±2ee separation. (16) 

We shall extend the investigation of this model over a 
broader range of parameters than suggested by its initial 
physical motivation. 

Note that, in the continuum limit, theories with high 
order derivative terms will generally give rise to 

v{k) = P(fc2) (17) 

where P is some polynomial. Although Vz{k) and its likes 
are artificial on the lattice, their continuum limit is quite 
generic. Later on we will show that if P(fc^) attains its 
global minima at finite \k\, then thermal instabilities can 
incur an extremely low value of Tc- 



In an "Ising" system S{x) = ±1 everywhere on the 
lattice. Stated alternatively, the scalar (n = 1) spins 
satisfy a normalization constraints 

{S\x) = 1} (18) 

at all A'' lattice sites x. Henceforth, we will adopt the 
latter point of view. 

Let us define the manifold M spanned by the set of 
minimizing wave-vectors q 

v{q€ M) = ia:m{v{k)}. (19) 
k 

If the local normalization constraints are swept aside 
then it is clear that the ground states are superpositions 
of sinusoidal waves with wave- vectors qGM. One would 
expect this to be true, in spirit, also in the highly con- 
strained Ising case, if v{k) is sharply dipped at its global 
minima. "Digitizing" a particular plane wave 



^^^^ IV. ISING GROUND STATES 
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S{x) = sign{cos{qi ■ x)) 



(20) 



and comparing it with the exact (numerical) ground 
state, one finds encouraging agreement in certain cases. 
For instance, this gives reasonable accord when H — 

jjMott 

This Hamiltonian (with some twists) was investigated 
in H] on a square {d = 2) lattice. 

Note that in the continuum limit (i.e. if the lattice is 
thrown away) we might naively anticipate a huge ground 
state degeneracy- a "digitized plane wave" for each wave 
vector flying on the (d— 1) dimensional manifold {Mq : 

= Q} This large degeneracy might give rise to a loss 
of stability against thermal fluctuations. 

It is found that striped phases (i.e. "digitized plane 
waves") were found in virtually all of the parameter 
range. Only for a very small range of parameters were 
more complicated periodic structures found. 

An intuitive feeling can be gained by considering a one 
dimensional pattern such as 

+ + -- + + -- + + --... (21) 



This pattern is a pure mode 

I— TX TT 

Speriod=Ax) = V2cos[-a; - -J. 
A double checkerboard pattern such as 



(22) 



+ 



+ - 



- + H 

+ H 

h + 

h + 



(23) 



extending in all directions in the plane is thus trivially 
given by 

.TT TT TT TT 

b\x) — 2cos — COS — a;2 = 

cos[-(xi + 2^2) ~ 2 ^ *^°^[2 ~ ^^^^ 

Such a 4 X 4 X 4 periodic pattern in three dimensions 
would include the eight modes |(±7r, ±7r. ± tt). This 
trivial example serves to illustrate an simple point. If one 
has a periodic building block of dimensions pi x p2 x Ps, 
then 



S{x) = Sj,\xx)Sj,^{xi)Sp^(xz). 



(25) 



If a configuration Sy{x) contains the modes {fc™} 
with amplitudes {S'p(fem)}, then Fourier transform- 
ing the periodic configuration S(x) one will find the 
modes (ibA:™^ , ±fc™^, zbA;™^) appearing with a weight ~ 
|5'pj(fci) X Sp^ik-i) X 5'p3(fc3)p. For high values of the pe- 
riods the weight gets scattered over a large set of wave- 
vectors. If i'(fc) has sharp minima, such states will not 
be favored. The system will prefer to generate patterns 



s.t. in all directions i albeit one Pi = 1 (or perhaps 2). 
For a, pi X p2 X P3 repetitive pattern, the discrete Fourier 
Transform will be nonzero for only II^^j^ pi values of k. 
This trivial observation suggests the phase diagram ob- 
tained by U. Low et al. jsj in the two dimensional case. 
The intuition is obvious. We have derived rigorously, 
the ground states in only several regions of its param- 
eter space (those corresponding to ordering with half a 
reciprocal lattice vector), and on a few special surfaces 
(corresponding to ordering with a quarter of a recipro- 
cal lattice vector) . In all of these cases the Ising states 
may be expressed as superpositions of the lowest energy 
modes exp[iqm ■ x]. Lately, a beautiful extension was car- 
ried out by 0. 

We now ask whether commensurate lock-in is to be 
expected. The energy of the Ising "digitized plane wave" 
on an L X L X L lattice where q — (gi,0, 0) with qi = 
27r/m, with even m, reads 

1 



E 



k 



i=i,3,...,m-i I exp[27rij/TO] - 1|2 
m 



E 

J — l,3,...,m — 1 



sin (TTj/m) 



(26) 



The lowest energy state amongst all states of the form 
considered is a possible candidate for the ground state. 

For the particular model long-range introduced above, 
it seems that for small values of Q, it might be worthwhile 
to have an incommensurate phase. This is, in a sense, 
obvious- all low energy modes are of very small wave- 
number and hence not of low commensurability. The 
energy 



16v{k = {2n+l)q) 



(2n-hl)27r2 



(27) 



For q ~ Qi/^'^+i) ^ 1, the higher harmonics k — 
{2n + l)q, do not entail high energies. For large val- 
ues of Q, g ~ 0(1), and v{k = (2n -I- l)q) can be very 
large if [(2n -|- l)q\ approaches a reciprocal lattice vector 
K. Under these circumstances it will pay off to have a 
commensurate structure; for a ui x U2 x ... x repeti- 
tive block only the modes k = 2Tr(^, ^, ^) wiU be 

populated (i.e. have a non-vanishing |5(fc)p)- the fer- 
romagnetic point [a reciprocal lattice point] will not be 
approached arbitrarily close- if that is not true weight 
will be smeared over energetic modes. Generically, we 
will not be expect commensurate lock-in in limg — > 
for any theory with a frustrating long range interaction. 
Although we have considered only striped phases (which 
have previously argued are the only ones generically ex- 
pected) , it is clear that this argument may be reproduced 
for more exotic configurations. 
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For finite range interactions it is easy to prove, by 
covering the system with large maximally overlapping 
blocks, that there will be a sliver about ^=0, for which 
we will find the ferromagnetic ground state. 

A polynomial in A(fc) will have its minima at A(/c) = 
const, i.e. on a (d-1) dimensional hypersurface(s) in k- 
space or at the (anti) ferromagnetic point. 

The kernel Vz {k) = zA^ — A has its minima {z > 0) at 

geM, : A(g) = min{-L,4d} (28) 

ZZ 

For z > -^i Mz is {d — 1) dimensional. 

We may divide the lattice into all maximally over- 
lapping 5 X 5 X ... X 5 hyper-cubes centered about each 
site of the lattice. 

Energy = ^ ^ ^ e{hypercube) (29) 

hypercubes 

and evaluate the energies e of all 5 x 5 x ... x 5 Ising 
configurations. Of all 2^ configurations the Neel state 
will have the lowest energy for a sliver about z = 
Analogously for z > Ztop » 1, by explicit evaluation, 
the ground state will be ferromagnetic. Contour argu- 
ments can be employed and a finite lower bound on 
generated. 

In this system on the square lattice with periodic 
boundary conditions along all diagonals e± : defined by 
xi^ X2 = const Ising ground states for 2 = | can be 
synthesized. All minimizing modes lie on 

M,^i : |(7i±g2| (30) 

By prescribing an arbitrary spin configuration along X- 
and fixing S{x_,x+) = S{x_,Q){-lY+: 

S{k) = ^5'(a;_,0)exp(ifc_a;_)^(-l)^+exp(ifc+a;+) 

X- X-l- 

(31) 

vanishes for = |fci + ^2! 7^ tt- Similarly by tak- 

ing the transpose of these configurations we can generate 
patterns having S{k) = unless \k-\ = \ki — k2\ = n. 
The ground state degeneracy is bounded from below by, 
the number of independent spin configurations that can 
be fashioned along x+ or x-, (2^+^ — 2) where L is the 
length of the system along the x± axis. The number of q 
values, commensurate with the diagonal periodic bound- 
ary conditions, lying on M^^x is {4L — 2). 

[Similarly for d > 2, one can set {d — 2) of the q com- 
ponents to zero. There are d{d — l)/2 cross-sections of 
the d-dimensional M^^i, all looking like the the two- 
dimensional M just discussed (i.e. \qi ± q2\ = n ). The 
real-space ground state degeneracy is bounded from be- 
low by d{d — l)[2^ — 1] (along the {d — 2) zero-mode 



A) FLIP SPINS DIAGONALLY 




FIG. 2. Simple ground states for «z(fc) with z = 1/8 (and 
for VQ{k) when Q — 16): A) Their construction- along the ar- 
rowed line we arbitrarily prescribe spins. For each move along 
the dotted diagonal lines we flip the spins. B) In this chain 
of diagonal super-spins there is no stiffness against flipping. 



directions the ground state spin configurations display 
no flip).] 

If we regard each diagonal row of spins as a "super- 
spin" then we will see that flipping any "super-spin" en- 
tails no energy cost. This is reminiscent to a nearest 
neighbor Ising chain where the energy cost for flipping 
a spin is dwarfed by comparison to the (logarithmically) 
extensive entropy. We might expect that here, too, or- 
dering might be somewhat inhibited. 

In two dimensions 

lim Mz:q'^ = — (32) 

z^oo 2z 

the "average" number of allowed q G Mz values <C 
0{L) (and similarly for the onset lim^^j^+ Mz : (^ — 

(±7r,±7r))^ = (g^ — ^)). For a hypercubic lattice of 
size Li X L2 X ... X Lg in d > 2 many discrete recipro- 
cal points will give rise to the same value of A(fc) ~ P. 
The proof is trivial: if all L; = L, then the number of 
possible k'^ values is bounded by dL"^, whereas there are 
L"^ fc— values. 

Therefore, on "average" , the number of k points lying 
on Mz, or more precisely lying lying the closest to Mz, 
s.t. |A(^) — ^1 is min, is, at least, 0{L'^~'^). 

[Of these, f^' wave-vectors, with rij (and z) denot- 
ing the number of identical components (and the number 
of zero components) of a certain fci nearest to Mz, are 
related to ki by symmetry.] 
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As we have stated previously, in the continuum limit 
(q 0) any short range kernel (including this one) will 
have a uniform (ferromagnetic) ground state. However 
the impossibility of constructing ground states that con- 
tain only "good" Fourier modes S{k S M) when M 
shrinks to a curved surface enclosing the origin is more 
general and will proved in the next section. 

For this short range model, even for z 7^ ^ a huge 
ground state degeneracy is expected. A "plane wave" 
might correspond to each wave-vector q (or commensu- 
rate wave- vectors nearby) lying on the (d—l) dimensional 
manifold M. 

As we shall prove later on, even in high dimensions, 
and even if the interactions are long ranged, in the con- 
tinuum limit it will not be possible to construct Ising 
states in which S{k ^ M) = unless the minimizing 
manifold M contains flat non-curved segments (or more 
generally intersects a plane at many points). 

V. A UNIVERSAL PEIERLS BOUND 

If a real hermitian kernel v{k) attains its minima in 
only a finite number of commensurate reciprocal lattice 
points {(ji}, then a Peierls bound can, in some instances, 
be proven for an infinite range model: When possible this 
is suggestive of a finite Tc. 

For instance, the bound for a (lattice) Coulomb gas 
(with the kernel solving the discrete Laplace equation on 
the lattice) is trivially generated. 

Ve{k) - Ve{q) = e/A(fc) - e/A{q = tt, 7r,7r) 

>-A(A(fc)-A(g = 7r,^,7r))), (33) 

with A = 16(P. The right hand side is the kernel of 
an antiferromagnet. Both system share the same ground 
states. For a given configuration the energy penalty for 
the Coulomb gas 

AEe = l/{2N)Y,[ye{k) - veimsm'' (34) 

k 

is bounded from below by the corresponding penalty in 
an antiferromagnet of strength A. In d = 2 the contour 
penalty of the antiferromagnet is 2A|r|, (|r| = length 
of the contour F). A similar trick may frequently be 
employed when the minimizing wave- vectors attain other 
commensurate values. It relies on comparison to a short 
range kernel for which a Peierls bound is trivial. 

VI. ISING WEISS MEAN-FIELD THEORY 

Assuming that for T < Tc- {S{x)) ~ s sign[cos{q ■ x)], 
then 



Y^ismvi^s = 0, y) - ^ E(^) (^)"(-^) 

y k 
= ' E (2n+ iw [^(^ = + m) + v{k = -{2n + l)q)] (35) 

n=0 ^ ' 

where we have assumed Qi = ^ with ^ 1 for all 
d components in replacing a discrete Fourier transform 
sum by an integral. 

s^{S{x^Q)) = -teinh[pY,V{x = 0,y){S{y))] (36) 

y 

yields 

n=0 ^ ' 

The latter is an upper bound on Tc as the non-trivial 
solution [with a non-zero magnetization s] should be self- 
consistent at all sites (not only at a; = 0) 

Summary and Outlook 

We have argued that for translation invariant two- 
spin kernels, the ground states are generically uniform, 
Neel ordered, or thin stripes (all periodic configurations 
in which the net weight |5'(A:)P is not smeared over 
many non-identical modes.) If the Fourier transformed 
kernel v{k) is sharply peaked about its global minima 
(at {q\) and if all other local minima have a much higher 
"energy" v{k) then the modulation length is 0(|(7|^^). 
All this might seem trivial/naive and incorrect. (It is 
also possible to get a feel for this via the examination 
of the local magnetization {S{x))) at low temperatures 
as computed within the dual unconstrained Hubbard 
Stratonovich Hamiltonian.) 

To emphasize that this is only generic but non- 
universal, we have constructed non-periodic ground 
states of a certain finite ranged Hamiltonian. 

As we have seen, frustrating interactions can give rise 
to massive degeneracy. These could, potentially, give rise 
to a small value of Tc. 

In the continuum limit [\q\ —t 0), the ground states in 
the presence of a frustrating long range interaction are 
expected to be non-commensurate; if short range frus- 
tration is present then the continuum limit ground state 
will be homogeneous (ferromagnetic). 

In this context, it might be worthwhile to point out 
that the triangular antiferromagnet was solved exactly 
by Wannier Q. This system has finite entropy at T = 0. 
There is no transition at finite temperature- this, again, is 
presumably related to this high ground state degeneracy. 

VII. 0{N > 2) GROUND STATES 

By an 0{n) system we mean that the spins {S{x)} have 
n components and are all normalized to unity- S'^{x) = 1. 
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Our previous ansatz S{x) = sign(cos{qi ■ x)) is readily 
fortified in the 0(n > 2) scenario: here there is no need 
to "digitize"- in the spiral 



Si {x) = cos(g • x) 
S2{x) = sm{q ■ x) 
Si>2{x) = 



(38) 



the only non-zero Fourier components arc S{q), S{—q). 
In plain terms, this state can be constructed with the 
minimizing wave-vectors only. 

It follows that any ground state g must be of the form 



(39) 



Let us turn to the normalization of the spins at all 
sites. 

n ^ n 

l(f)^5^5f(f) = -5]5]arar' 

(C0s(^r + 4>f)cOs[(qm + q-m') ' x] 

- sin((/)™ + 0„/) sin[(g;„ + (f„/) • x] 
+ cos{4>^ - 4>f )cos[{qm -qm')- x] 
- sm{<Pr ~ €'')M{q^ - Qm') ■ x]) (40) 

If l(x) = 1 is to hold identically for all sites x, then all 
non-zero Fourier components must vanish. 
For the {cos(A • x)} Fourier components: 

n 
n 

+ E E«rar'cos(C-C')] (41) 

and a similar relation is to be satisfied by the {sin(^ • 
x)} components: 

n 

0=[ E E<"r'sin(C + C') 



3^1 =A » 



= 1 



+ J2 E«>"'^i'^('^"-'^r')] (42) 

(43) 



a^cosC^^r, a™sinC = <- 
The {cos(^ • x)} and {sin(A • x)} conditions read 

n 

0= E E[«'-<x"'] 

n 

+ E 

n 

0= E E[<«r'+«>r'] 



+ Yl E[<- - 

In the case of two pairs of wave-vectors ±qi and ±92, 
both not equal half a recipro- 
cal lattice vector: (0, 0, .., ), (tt, 0, 0), (0, TT, 0, 0),..., 
(tt, 7r,0, ...,0), (tt, TT, ...,7r), the vector A (up to an ir- 
relevant sign) may attain four non-zero values: A = 
2qi,2q2,qi^± 9*2 • 

When A = qi + q2, the conditions are 

n 

71 

= E[^N? + «k']- (44) 

i=l 

When A = qi — q2, these conditions read 

n 

o = EK'«'+^'«i] 

i=l 
n 

o = YMy^-^^] (45) 

i=l 

Fori* =24 (a = 1,2): 

n 

= E[« - «] 

i=l 

n 

= 2^« (46) 

i=l 

Define 



f/« = (u;ii, <,-,<) 
F« = «,...,<) 

The previous conditions imply that 



(47) 



(48) 



The four vectors {U^ jV"^} are all mutually or- 
thogonal. The number of spin components n > 4. 
Two additional demands that follow are 

2 

^[F" • -|-f7" • i/"] = 2^\7" • = 2. (49) 

a=l a 

The last equation is the normalization condition- the 
statement that the coefficient of cos(A • x), when A = 0, 
is equal to 1. 
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For the case of a single pair of wave- vectors, ±gi , A ^ 
2qi , and the sole conditions are encapsulated in the last 
of equations and in equation (^). 

A moment's reflection reveals that this only allows for 
a spiral in the plane defined by and . 

When n < 4 there are no configurations which satisfy 
S'^{x) = 1 identically for all sites x [excusing those having 
2{qi + qj) = Ais equal to a reciprocal lattice vector] that 
are a superposition of exactly two modes. 

For instance, a (double checkerboard state along the 
i — 1 axis) (E" (a spiral in the 23 plane) has pairs {i,j) 
for which A = 2{qi + qj) is a reciprocal lattice vector. 

As the number of minimizing modes {qm} increases, 
some of the conditions may degenerate into one, e.g. if 
{qi + Q2) — {(l3 — (I2) (i-e. the modes are colhnear). This 
degeneracy is the a second route that might allow for 
Ising configurations which are superpositions of several 
"good" minimum energy modes exp(i(f • x). 

The highly degenerate Ising ground states that we have 
constructed previously can fall under either one of these 
categories. 

If neither one of these situations occurs, Ising states 
cannot be superpositions of several minimum energy 
modes: we will be left with too many equations of con- 
straints with too few degrees of freedom. 

For three pairs of minimizing modes, none of which is 
half a reciprocal lattice vector, {±qm}m=i with 

qnj±qt^qr±qs^ 2qw (50) 

for all w ^ t, and r 7^ s, conditions similar to those 
that previously written for yl = gi ± 52 , now hold for all 
[q^ ± qt). 

Ua ■ Up — Sa,l3 

Uo.-% = ^ (51) 

The relation U a • Va — ^ [ol = ^ va the last eqn above) is 
enforced by setting A — 2qa- Thus, when exactly three 
pairs of minimizing wave- vectors satisfying the equation 
are present, the vectors {Ua,Va\ define a 6-dimensional 
space, and hence n > 6. For p pairs of minimizing wave- 
vectors, n must be at least 2p— dimensional. This bound 
is saturated when S is (a spiral state in the 12 — plane) 
(g) (a spiral in the 34 — plane) (a spiral in the 

2p— l,2p plane), i.e. 

[ax cos((fi • X + (t)i),ai sin((f • x -f- 0i), ... 
, Op cos(qp • x + (j)p),ap sin(qp • x -|- (f)p)) (52) 

with '^a — ^- When wave- vectors with q^ ± qt = 

qr ± qs or (lui ± qt — 2qr are present, pairs of conditions 
degenerate into single linear combinations. 

So far we have assumed that for all i and j, A = 2{qi + 
qj) is not a reciprocal lattice vector, s.t. sin(A • x) is not 
identically zero at all x G Z'^. 



We term the such a p = 2 configuration a bi-spiral. 
It is simple to see by counting the number of degrees of 
freedom for n = 4, that the bi-spirals overwhelm states 
having only one mode ±gi. This is a simple instance 
of a general trend: High p states are statistically pre- 
ferred. Moreover, as we shall see later, they are more 
stable against thermal fluctuations. 

Summary and Outlook 

We have outlined a way to determine all 0{n > 2) 
ground states for a given kernel V{x,y) = V{x — y). 

Whenever n > 2, any ground state configuration 
can be decomposed into Fourier components, S^{x) = 
Sl^i {cos[qi ■ x]+bi sin[(fj • f]| 

qi are chosen from the set of wave vectors which mini- 
mize v{k). 

\M\ is the number of minimizing modes (the "mea- 
sure" of the modes on the minimizing surface M . 

So long as these wave- vectors qi which minimize v{k) 
are "non-degenerate" , in the sense that the sum of any 
pair of wave vectors, qi ± qj is not equal to the sum of 
any other pair of wave vectors, and "incommensurate" 
in the sense that for all i and j, 2{qi + qj) is not equal 
to a reciprocal lattice vector, then the that the condition 
[S'f (f)]^ = 1 can be satisfied only if \M\ < n/2. (In our 
toy model of the doped Mott insulator, these conditions 
are always satisfied for Q < 4.) Thus, for n < 3 only sim- 
ple spiral — 1) ground-states are permitted, while 
for n = 4, a double spiral saturates the bound. Thus, 
generically, for 2 < n < 4 all ground states will be spirals 
containing only one mode. 

The reader should bear in mind that in the usual short 
range ferromagnetic case, the ground states are globally 
SO{n) symmetric and are labeled by only (n—1) contin- 
uous parameters. 

Here, for each minimizing mode there are (2n — 3) con- 
tinuous internal degrees of freedom labeling all possible 
spiral ground states. For n > 2 this guarantees a much 
higher degeneracy than that of the usual ferromagnetic 
ground state. 

If there are many minimizing modes (e.g. if the mini- 
mizing manifold AI were endowed with SO{d — 1) sym- 
metry) then the ground state degeneracy is even larger! 

When n > 4, there are (generically) even many more 
ground states (poly-spirals). These poly-spiral states 
have a degeneracies larger than those of simple spiral. 
Their degeneracy 

g^pi2n-2p-l)\M\P, (53) 

where |M| is the number of minimizing modes. 

To capitulate: we have just proved that if frustrat- 
ing interactions cause the ground states to be modulated 
then the associated ground state degeneracy (for n > 2) 
is much larger by comparison to the usual ferromagnetic 
ground states. 
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VIII. SPIN STIFFNESS 

When Q > the minimizing modes He of Vcontinuum {k) 
lie on the surface of a sphere {Mq : cp = \/Q}- 

As Q ^ 0, this surface Mq shrinks and shrinks yet is 
still a (d—l) dimensional surface of a sphere. When Q = 
0, the minimizing manifold evaporates into a single point 
^ = 0. This sudden change in the dimensionality has 
profound consequences. As we shall see shortly, it lends 
itself to suggest (quite strongly) that order is inhibited 
for a Heisenberg (n = 3) realization of our model. 

Before doing so, let us indeed convince ourselves, on 
an intuitive level, that the large degeneracy in k— space 
brought about by the frustration gives rise to a reduced 
spin stiffness. 



A. Longitudinal 



Let us assume twisted boundary conditions 

S(x) = cos(— a; + qx)ei + sm(— a; + qx)e2 + oS 

= ei cos(fc • ^) + 62 sin(^ • x) + 55 (54) 



t;(fc) — vi^q) = 



167r-^ 



Ignoring SS contributions: 



AH = 



AH ^0 



(61) 



(62) 



(63) 



as L — > oo in d = 3 [a complete loss of stiffness against 
transverse fluctuations]. For a one dimensional chain of 
nearest neighbor XY spins: 



(64) 



For our frustrated three-dimensional system, the trans- 
verse fluctuations obey 



(65) 



with k = {^+q}ei. The energy cost of this state relative 
to the ground state is 

AH[{Six)}] = ^ ^[z;(fc') - v{q-)] \S{k')\' (55) 



Ignoring 6S contributions: 

AH = -Hk) - vm = ^[(^ + if - - 



For the usual nearest neighbor ferromagnetic XY model: 
S{x) = cos[^]ei + sin[^]e2. (57) 

Here 



v{k) = 2^(1 - cos ki) 



1=1 



AHxY = [l-cos{^)]N-^^ 



(58) 



(exactly the same). 



In principle, one may envision other impositions of 
twisted boundary conditions [say in the 3-4 plane]: 

S{x) = [cos{qx)ei + sin(ga;)e2]\/l — 



(56) 



+ [cos(^)e3 + sin(^)e4]e + 55 (66) 



AH=^[v{k='^ey)^v{g)]N 



2^2 



= ^e^[v{k = 0) - v{q)] = ^ 



(67) 



where in the last line we have taken the limit L ^ oo. 
The system exhibits infinite spin stiffness to these sorts 
of fluctuations. 



D. 



To summarize, the system responds to fluctuations 
with an effective kernel 



B. Transverse 

— * ^TTX 2tTX —* 

S{x) = cos[^— + qy]ei + sin[— — + qy]e2 + 5S (59) 

Ij Ij 



k = — ei +qe2. 



(60) 



(68) 



where S± and denote the the transverse and longitu- 
dinal fluctuations. 

For a nearest neighbor one dimensional system embed- 
ded in d-dimensions, v{k) — v{q) ~ ^W^'^i (^x = 0), and 
thus the fluctuations are even larger than in our case. 
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IX. THERMAL FLUCTUATIONS OF AN XY 
MODEL 

Let me begin by treating the "soft-spin" version of the 
XY model, in which we include the non-linear interaction 



H. 



soft 



Ho 



X 

= Ho + H, 



(69) 



with u > and forget about the normalization conditions 
\S{x)\ = 1 at all lattice sites x. 

(The normalized "hard-spin" version can be viewed as 
the u — > cxD limit of the soft-spin model.) 

As we have seen previously, the only generic ground 
states (for both hard- and soft-spin models) when the 
spins S{x) have two (and also three) components are spi- 
rals 

ciQround — state / ^\ r -> -'\ r^qround — state t ^\ • ^\ 

Df (a;) = cos(g • a;); d| [x) — sva[q ■ x) . 

(70) 

We will expand Hsoft about these ground states, keep- 
ing only the lowest order (quadratic) terms in the fluc- 
tuations 5S. The quadratic term in {5Si{k)} stemming 
from Hsoft is the bihnear j^{5S)^ M{8S) where 

{5S)+ = (<55i(-fci),J52(-fci),<55i(-fc2)<552(-fc2), 

5Si{-k^), 6S2{-h), SSi{-kM),SS2{-kN)) (71) 

and the matrix M reads 



/ 4 . 


. 1 


i 


. . . N 


4 . 


i 


-1 . 




1 -i . 


. 4 


. 


1 i 


-1 . 


. 


4 . 


. i -1 




. 1 


— i 


. 4 


V ! 


—i 


-1 . 


.0 4 / 



A few off-diagonal sub-matrices have been omitted. 
[There are two along each horizontal row (by periodic 
b.c.).] The sub-matrices are (2 x 2) matrices in the inter- 
nal spin indices. The off diagonal blocks are separated 
from the diagonal ones by wave- vectors (±2q). 

Note that (fc|M|fc') = M{k - k'). Making a unitary 
(symmetric Fourier) transformation to the real space ba- 
sis: 



(72) 



the matrix M becomes block diagonal 



{x\M\x')=M{x)5s,s'- 
Diagonalizing in the internal spin basis: 
A± =6,2. 



(73) 



(74) 



These eigenvalues may be regarded, in the usual fer- 
romagnetic case {q = 0) as a two step (state) potential 
barrier separating the two polarizations. I.e., the nor- 
malization constraint of the XY spins (embodied in M) 
gives rise to an effective binding interaction. As we shall 
later see, when the number of spin components n is odd, 
one spin component will remain unpaired. Hi literally 
"couple" s the spin polarizations. 

Employing Eqn.([74|), we note that the corrected fluc- 
tuation spectrum {ipm]m=i (to quadratic order) satisfies 
a Dirac like equation 



[U+ v{~id^) U + 2u 



2 
6 



Em U+\^Jm{x)), 



Here 



U 



sin(2q-x) l 
2cos(g-:r) 2 

— 1— cos(2g-^) 1— cos(2g-:E) 

2 cos{q-x) 2 sin(q-x) 



(75) 



(76) 



Alternatively, expanding in the fluctuations 6S{k). 
leads to bilinear {SS)'^H{SS) where 



v{k) + 8u 
v(k)+8u 



(77) 



along the diagonal, and 



2u 2iu 

2iu —2u 



2u 



-2iu 



'2iu —2m 



(78) 



off the diagonal. 

^fc,fc±2g (a3±«fTl) 



(79) 



exp(ij o-i)Hfc,fc±2g-exp(-i J ai) 



= 2u(T^ = Au [ 











r 






ol 



r = [v(k) + 8u] 



1 
1 



(80) 



If Q = 16, g = (7r,0,0) G Mq, and 2q = 0{mod 2tt). 
The fluctuation matrix is diagonal in the k basis, and 
the fluctuations are divergent at flnite temperatures. An 
identical situation occurs for q = (tt, 7r,0) € Mq=64- 
When Q = 4, g = (^, 0, 0) e Mq, and it is easy to show 
that determining the eigenvalue spectrum degenerates 
into a problem in two parameters (A2,fci) where A2 = 
2Y,^i^^{^~coski) {s.i. A{k + q) ^ A2 + 2coski, A{k) = 
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A2 — 2 cos ki ) , The fluctuation integral about the chosen 
ground state exhibits a (d — 2) dimensional minimizing 
manifold (parameterized, in our case, by (A2, fci) ). Note 
that, at higher order commensurabilities, the dimension- 
ality of the minimizing manifold is low. In fact, the in- 
teraction will no longer be diagonalized in fc— space. 

Till now, all that was stated, held for arbitrarily large 
u- our only error was neglecting 0{{SS)^) terms by com- 
parison to 0{{5S)^). Note that the main difficulty with 
the approach taken till now was the coupling between k 
and k ± 2q: i.e. k is coupled to A; ± 2q, while k + 2q 
is coupled to k + iq and k, and so on. Unless q is of 
low commensurability an exact solution to this problem 
is impossible. 

To make progress let us assume that u is very small. In 
this case the lowest eigenstates of the fluctuation matrix 
will contain only a superposition of the low lying k- states 
[i.e. those close to the {d — 1) dimensional M {Q > 0)]. 
li ki = q + S is close to M, then the only important 
modes in the sequence {Si{k = ki + 2nq)} are ki, and 
k2 = ki — 2q = —q + S. The sub-matrix in the relevant 
sector reads 

/'v{ki) + Su 4u \ 

v{ki) + 8u 

v{k2) + 8u 

V 4u v{k2) + 8uj 

The lowest eigenvalue reads 

Blow - ^[v{ki) + vfy] + Au 



v{ki) ~ v{k2)]^ + Giu^ 



(81) 



Equivalently, this can be determined from the direct 
computation of the determinant to 0{u^): to obtain 
0{u^) contributions we need to swerve off the diagonal 
twice. 



+ (4m)4 ^[w(4) + 8u]b(4 + 4?) + 8u] 
jl 

(The partition function is trivially 
Z = const [det7^]-i/2.) 

For any u, no matter how small, there exists a neigh- 
borhood of wave- vectors k near q such that \v{k) — v{k + 
2q) I ^ u and as before we may re-expand the character- 
istic equation for these low lying modes, solve a simple 
quadratic equation, expand in the components of (fc — (7) 
and obtain a simple dispersion relation. 



If (5 = (5_L + (5|| e||, with _L 



denoting directions 



orthogonal, and parallel to n _L Af , then 



where 



A 



II 



l (Pv{\k\) 
2 dk^ 



\k\=q] Aa 



4g2 



(83) 



(84) 



We have already found this form when examining spin 
stiffness of the "hard" XY model. (In both cases there is 
also an identical ^j^^i term which may be omitted.) Hav- 
ing reassuringly re-derived this dispersion from another 
vista, let us note this dispersion is akin to the fluctua- 
tion spectrum of the smectic liquid crystals ||l7|] which is 
well known to give rise to algebraic decay of correlations 
at low temperatures. In our case, one notes that q ^ 
and thus the correlations should have an oscillatory pref- 
actor. To be more precise, the correlator, in cylindrical 
coordinates, 



4^2 

G{x) 2- exp[-2?77 - riEi{' 



2x11 



)] X cos[qx\\] (85) 



detn 



where 77 — ^^^^Q^^^, ^ X "^here A is the ultra violet 



'Hi.ii'H2,i2---T^2N,i2N ■ (82) momentum cutoff, 7 is Euler's const., and 



Il,l2,-..,I2JV 



det n = n^^i [v{ki) + 8uf - {Auf ^ [v{kj) + 8u] 

j 

x[v{kj + 2q) + 8u] %^fe^,fe^+2«b(^.) + 

The fluctuation spectrum is trivially determined by re- 
placing v{k) by [v{k) — E] in det Ti and setting it to zero. 
To this order we re-derive Eiow ■ To higher order 

. . . + {Auf ) + N [^(^Ji + + 8"] 

[w(%) -I- 8u][w(% +2q} + 8u] 



i?i(z) = -7-ln2-^(-) 



^ n{n\) 



(86) 



is the exponential integral. For the uninitiated reader, 
we present this standard derivation in the appendix. 

The thermal fluctuations / d'^k/Eiow{k) diverge as 
[— In |e|] with a e a lower cutoff on |fc — g| (|]. Such a log- 
arithmic divergence is also encountered in ferromagnetic 
two dimensional 0{2) system if it were exposed to the 
same analysis. Indeed, both models share similar char- 
acteristics albeit having different physical dimensionality. 

summary and outlook 

We have found that if Hgoft is indeed soft (u ^ 1) 
then, within our continuum XY systems, only quasi long 
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range (algebraic) order will be observed at low tempera- 
tures in the presence of generic long range (e.g. VQ{k)) 
frustration These systems essentially display smectic 
phases. 



Some smectic trivia 

Smectic A ordering involves a displacement u of the 
smectic along the z-direction (the direction of the molec- 
ular axis). It is customary to denote this displacement 
by u{x). As we shortly see, the kernel for the smectic is 
identical to ours. 

When the wave-vector k is along the z-direction the 
displacement is longitudinal, and the energy is of the 
elastic form ii3fcj^|u(fc)p, where B is the compressibil- 
ity for the smectic layers. When k is normal to z, the 
displacement is transverse and the layer separation. No 
second order in /cj_ the displacement costs no energy. In 
this case the restoring force is associated with a director 
splay distortion. In this case, the elastic energy is density 
is ■ fiY with K the splay constant.. As n is normal 
to the layers, one has 5n = — V_lu and an elastic energy 



\Kki_\u{k)\'^. Thus the kernel 



v{k) 



Bkl 



(87) 



As seen in the form that we obtained for G{x), the pen- 
etration depth A = ^ KjB determines the decay of an 
undulation distortion (splay director distortion) imposed 
at the surface of the smectic. 



X. A GENERALIZED 
MERMIN-WAGNER-COLEMAN THEOREM 
(APPLYING TO ALL ANALYTIC (IN A SPACE) 
INTERACTIONS IN TWO DIMENSIONS) 

We will now slightly generalize the Mcrmin-Wagner- 
Coleman theorem Q , : 

All systems with translationally invariant two-spin in- 
teractions in two dimensions with a real analytic Fourier 
transformed kernel u(fc) show no spontaneous symmetry 
breaking. 

Our approach is the standard one. We will merely keep 
it more general instead of specializing to ferromagnetic 
order and or interactions of one special sort or another. 

The magnetic field 



h(x) — hcos{q- x)e„ 



applied by itself would cause the spins to take on their 
ground state values. 

If n = 2 the unique ground state {s^{x)) to which a low 
temperature system would collapse to under the influence 
of such a perturbation is 



Sf{x) =sm{q-x). (89) 
When n = 3 the ground state is not unique: 



S!<ni^) = nsm{q-x) 

n-l 
i=l 



(90) 



and a magnetic field may be applied along two directions, 
with all the ensuing steps trivially modified. 
With the magnetic field applied 

n 

^ = 9 E E Vix;-y)SMS^iy) - K{x)Sn{x). (91) 



Note that the knowledge of the ground state is not 
imperative in providing the forthcoming proof ||ll|] . 

The standard idea [|l^ that we are about to exploit is 
the rotational invariance of the measure. 



d^ji I Ilsd"S{x)d{S^{x)-l)e-^" (92) 



(93) 



The generators of rotation in the [af3] plane are 



Lal3 = Sa 



d 



S. 



d 







de 



d'^S sis^ ~ 1) 



f{Si, Sa cos9 + Sf3 sin 6*, ... 
, Sp cos9 — Sa sina, S'„). 



= / d" 5,5(5^ - l)Lapf{S). 



(94) 



(95) 



Now let us consider in particular the generators of rota- 
tion from the axis of the applied field to any other internal 
spin axis 



^ d d 

{Ls)i = Sn{x) - Si{x) 



dS,{x) 



dSn{x)- 



(96) 



In the up and coming _L will denote the projection per- 
pendicular to h (i.e. in the i < n subspace.) 
Let us define the following operators 



A{k) = ^exp[ifc • x]S±{x) 

X 

B{p) = Y,^Mw-ALsm). 



(97) 



By the Schwarz inequality. 

KE^r^oi' < (Ea*a.)* (5^5*5.). (98) 
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For any functional C: 

Ls{e-^"C) = e-f'"{L3{C) + CL(-/3H)}. 

= Jus d^'Six) 6{S^{x) - 1) L[e-^"C] 
{CB{p)) = (X^exp[jp. x] Ls{C)) 

X 

Lym] = 2 X 1/2 X /3^{[5„(y)5,(z) 

Z 

-Si{y)S„{z)]V{y~ z) ~ hn{x)Si)x)}. 

n n — 1 



i=l 



-h{x)[n-l)Sn{x)]) 
{B{p)* ■B{p))=(3Y.{{cosp-{x^y)-l) 

x,y 

n 

4=1 

V{x-y)}~{n~l)h{x){SrXx)) > 
Fourier expanding the interaction kernel 

t 

and substituting 

{S{x)-S{y))^^Y.^\S{u)\ 



we obtain that 



{B{p)* ■ B{p)) = ^Y. +P)+v{u-p) 



{A{kr-A{k))^Y.^s^{x)-SAy)) 



X exp[ifc • (a* — y)]. 



{A{kr ■ m) = {Y,LUS^{x))eMip- x]) 



(99) 
(100) 

(101) 



(102) 



(103) 



(104) 



(105) 



(106) 



2v{u)\{\S{u)\')-h,^{S^{~q)) (107) 



(108) 



where q = p—k and = {Sn{q))- Note that with our 
convention for the Fourier transformations, a macroscop- 
ically modulated with wave-vector q state would have 
m, = 0(iV). 

The Schwarz inequaUty reads 

2\m,f[pY.({\Siu)\^)ivip + u) 



(110) 



+v{p~ u)) - 2v{u)] + 2{n ~ l)\h\\mg\ 

<iV-i5](5x(f)-5i(y)). 



In the thermodynamic hmit: 



2(n-l)2 



d-^p 



{27r)nj {27rr 
{\S{u)f){v{p + u) + v{p- u) - 2v{u)) 

+2{n-l)\h\\mq\ 

<N^Slix)). 



However, 



{Slix)) < 1. 



(Ill) 



(112) 



In Eqn.(lll) p takes on the role of the deviation i5 
introduced in the earlier sections. Notice that in the low 
temperature limit. 



\p\>s d-^p 



(2^) 



(113) 



is non-negative 



Explicitly, as the integral J 

(as {B{p)* ■ B{p)) > the denominator in Eqn.(llO) is 
always positive for each individual value of p) , 

2(n- 1)2 



P\<S 



\mq\ 
d'^u 

{2nY L J {2nY 
{\S{u)f){v{p + u) + v{p- u) - 2v{u)) 

+2in-l)\h\\mg\ 

< N'iSlix)) 



(114) 



Taking S to be small we may bound 



p\<S jd 



d^p 



d<^u 



(n - l)m,- (109) 



{2tj:Y L J {2'kY 
{\S{u)\'^){v{p + u) + v{p- u) ~ 2v{u)) 
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+2{n - l)\h\\mq 



< 



\p\<s 



(fit 



{2t:Y L J {2ttY 
A,p'\s{\S{u)\^) 



+2(n - l)|/i||TOq| 



(115) 



with chosen to be the largest principal eigenvalue 
of the d X d matrix didj[v(u)], and Ai a constant. 
For analytic v{u), and for |p| < (5 where S is finite, 

iv{p + u) + v{p- u) - 2v{u)) < AiXsp"^ < Bip^ (116) 

for all u in the Brillouin Zone with positive constants Ai 
and Bi. 

In c? < 2, the integral 



(117) 



diverges making it possible to satisfy eqn.(lll) when the 
external magnetic field /i — > only if the magnetization 



■niq = 0. 



The finite temperature behavior of a quantum system 
is, in certain respects, similar to that of a classical system. 
One could directly tackle the quantum case by applying 
the Bogoliubov inequality 



^({A,At})*([ [C,i/],Ct])>|([C,^])p 



(118) 



with [ , ] and { , } the commutator and anticommu- 
tator respectively. Setting A = S^{k) and B = S'"(p) 



with q = p— k we will once again obtain equation(lll) 
with the classical spins replaced by their quantum coun- 
terparts. 



XI. 0(iV = 3) FLUCTUATIONS 

Let us now return to the more naive "soft-spin" 
0(3) models in order to witness an intriguing even-odd 
binding-unbinding effect that could have otherwise been 
missed. 

As we have proved, for an n = 3 system the generic 
ground states are simple spirals. If we rotate the helical 
ground-state to the 1 — 2 plane; the single quadratic term 
inSS,=3ix) is EsiSS^=3ix)Y. 



A 



j>3 



= 2 = A_ = A,i 



(119) 



Here, all that follows holds for arbitrarily large u- the 
only approximation that we are making is neglecting 
0{{SSY) terms by comparison to quadratic terms: i.e. 
assuming that dS{x) <^ 1. Unlike the above treatment 
of the XY spins, no small u is necessary in order to make 
headway on the Heisenberg problem. 

For n = 3 we find that the fiuctuation eigenstates of 
the are the products of an eigenstate of H within the 



plane of the spiral ground state and a fluctuation eigen- 
state of v{k) along the direction orthogonal to the spiral 
plane. Written formally, to quadratic order, the fluctua- 
tion eigenstates are 



\i^m)<E>\5S3{k)). 



(120) 



Fluctuations along the i — 3 axis are orthogonal (in a 
geometrical and formal sense) to the ground-state plane. 
This is expected as fluctuations in any hyper-plane per- 
pendicular to the [12] plane do not change, to lowest or- 
der, the norm of the spin. \SS3{k)) is literally the "odd" 
man out. As foretold, this is a general occurrence. When- 
ever the number of spin components is odd, one unpaired 
spin component is unaffected by the interaction enforcing 
the spin normalization constraint. 

Within the i = 3 subspace {k\M\k') = 26^ p and 
our previous analysis follows. The dispersion Elr^ — 
[v{k) + 2 u Ai^a] does not have a higher minimum than 
with [|l6|. 

Both have the same A value and the 6Si=3{x) fiuctu- 
ation is minimized at wave-vectors £ e Mq s.t. v{£) — 
v{q) — imiij:v{k). As \ipm) is a normalized superposition 
of \SS{k)) modes (the latter spin vectors being in the 1 — 2 
plane) and v{k) is diagonal, and attains its minimum at 
£& Mo: 



unn{E„,}>Egl^^. 
Here, we invoked the trivial inequality 
iRm{^/j\[Ho + Hi]\tP) > min((?!)|i7o|' 



(121) 



mm(e|i/i|e). (122) 



Thus, there exist Goldstone modes corresponding to 
SSi^3 fiuctuations, and one must adjust additive con- 
stants s.t. mmj:{E^^^} = 0. 

The resulting fluctuation integral reads 



1 



d^k 



(27r)3 v{k)-v{q) 



(123) 



(where we noted translational invariance 

(A5,(fc)A5,(fc')) =%+£,,o(l^'5»(^)l') and employed 
equipartition). As pointed out earlier, when Q > in 
1^ continuum (k) the minimizing manifold is {d — 1) dimen- 
sional. The fluctuation integral receives divergent contri- 
butions from the low energy modes nearby. By quadratic 
expansion about the minimum along n _L Mq, a diver- 
gent one-dimensional integral for the bounded {{AS{x — 
0))^) signals that are quadratic fluctuation analysis cal- 
culation is inconsistent. We are led to the conclusion 
that higher order constraining terms are imperative: We 
cannot throw away cubic and quartic spin fluctuation 
terms {6S^'^{x)) relative to the quadratic {6S'^{x)) terms 
(notwithstanding the fact that all of these terms appear 
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with 0{u) prefactors irrespective of how large u is). The 
spin fluctuations SSi^3{x) are of order unity at all finite 
temperatures and Tc{Q > 0) ~ 0. 

Note that this is "almost a theorem" . Here we do not 
demand that u be small (only SS). This is an important 
point. \ipm) is an eigenstate for arbitrarily large u. To 
quadratic order in the fluctuations, the minimum belongs 
to ISSs) (or is degenerate with it). 

The divergent fluctuations here signal that 0(55*^) — 
0{SS'^). Thus assuming that SS < y/jju (with J = 1 
the exchange constant) we reach a paradox. Thus, if the 
integral in Eqn.( p^3|) div erges for an 0(3) model then at 
all finite T: 6S > Jj/u. 

When Q — the minimizing manifold shrinks to a 
point- the number of nearby low energy modes is small 
and our fiuctuation integral converges in d > 2. This is 
in accord with the well known finite temperature phase 
transition of the nearest neighbor Heisenberg ferromag- 
net: Tc{Q = 0) = 0(1) (or when dimensions are fully 
restored- it is of the order of the exchange constant). 

Notice that a discontinuity in Tc occurs as Mq 
Mq^o = iq = 0). 

We anticipate that small lattice corrections (A 7^ 0) in 
VQ{k) to yield insignificant modifications to Tc{Q): one 
way to intuit this is to estimate Tc by the temperature at 
which the fluctuations, as computed within the quadratic 
Hamiltonian {5S^[x = 0)) = 0(1). 

Summary and Outlook 

We have argued that in (essentially hard spin) Heisen- 
berg realizations of our models no long range order is 
possible in the continuum limit: Tf.{Q > 0) = 0. More 
generally we claim that if the integral 



v{k) — v{q) 



(124) 



diverges then no long range order is possible at finite 
temperature. If lattice effects are mild then Tc{Q > 0) 
is expected to be small. In the unfrustrated {Q — 0) 
Heisenberg ferromagnet in d — 3: Tc — 0{1). 
Fusing these facts, a discontinuity in T^. 



STc = [7^c(0) - TciQ)] 



(125) 



is seen to exist. Tc{Q = 0) is an avoided critical point 
(see fig. 1 in the introduction)). 



XII. 0{N > 4) FLUCTUATIONS 

The fluctuation analysis of any 0{n > 2) system about 
a spiral ground state is qualitatively similar to that of the 
Heisenberg system. 

For a vanishing lower cutoff e on ||fc| — g|, the fluctua- 
tions of an n — component spin (d ~ 3) about a helical 
ground state are given by 



knT 



(126) 



More generally, in d > 2 dimensions, the leading order 
infrared contribution reads 



,d-3Ql/4 



(n~2) g('^-i)/4 

2d T,d/2 g T{d/2) ~ 2^+1 n(d-i)/2 (d-3) r(^ 



(127) 



As noted previously, poly-spiral states will tend to 
dominate at large n. 

The reader can convince him/herself that for even n 
with a, p < n/2 poly-spiral ground state and for all odd 
n the fluctuations will give rise to a leading order di- 
vergence. The reasoning is simple: the poly-spiral states 
extend along an even number of axis. If n is odd then 
there will be at least one internal spin direction i along 
which Sf = and our analysis of the Heisenberg model 
can be reproduced. 

The lowest eigen-energy associated with the fluctua- 
tions \ipm) in the {2p) dimensional space spanned by the 
ground state is higher than the lowest eigen-energy for 
fluctuations along an orthogonal direction. 

{iplHsoftW > 0. If |(55'| < 1, this implies that 
the quadratic term in SS{x) stemming from Hi is non- 
negative deflnite. 

For Sf{x) = 0, this quadratic term in 5Si{x) is zero. 

The eigenvalue Amin A_ = 2 corresponds to the zero 
contribution in 0{5S'^{x)) from Hi. 

And once again 



m.\Yi{E„,} > Ei, . 



(128) 



from the trivial inequality 

min{ip\[HQ + H^ofi 



> 



min(0|i7o|(/)) + imn{S,\Hsoft\C)- 

<t> C 



(129) 



The fluctuations of even component spin about a p = 
n/2 poly-spiral ground state are more complicated. 

Once again coupling between different modes occurs. 
In this case they arc more numerous. 

For rt = 4, the fluctuation energy, to quadratic order, 
about a bi-spiral reads 



5H 

^{Efc, ,-,J(55i(fci)<553(fc2) + <5^i(A:2)<553(fci)) 

^ gi+g2+fci+'£2,0 gi+92-fci— K2.0 fei+fc2+i32-gi,0 
%i+fe2+rfi-g2,o] + 

+ {6S2{k2)SSi{ki) + 5S2{k)5Si{k2)) 

I- Ki+K2+gi— 92.0 Ki+K2+g2— 91.0 fel +A:2 +91 +92 ,0 

%i + fe2-gi-92,o] 



i{5Si{k2)5Si{ki) -f 5Si{ki)5Si{k2)) 

X [%i +£2+91+92,0 ~ '^fci+fe2-9l-<?2,0 + %i+fe2+92-9l,0 
'^fei+fc2-92+9l.o] 

i{5S2{k2)5S3{ki) + SS2{ki)5S3{k2)) 
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^ ['^fel+fe2+9l+92,0 %i+fc2-9l-92,0 + %i+fc2+gi-92,0 



%i+fe2+<?2-«l,o]]j' 



2ua 



-E£,,£,{<55i(^i)<55i(fc2)[<5,^,+£„o + 5(%.+fc2+2,i,o + 



'^fcl+/C2-2qi,0>' 



%i+iS2-2gi,o)]} 



2ua2 



)] 



(554(A;i)(5S'4(A:2)[(5g^_^^^_o 



'^fel+fe2-2<f2,o)]j' 



2V"fei+/c2+2g2,0 



^E£„&[^^i(fc2)^52(A;i) + 

- ^ ^ ¥E£„fe[^^3(^2)'554(fcl) + 

(5S3(fcl)554(fc2)](%^+fc,+2,-2,0 - h,+k2-2q2fi)- 

As before, wc may obtain an equation for the eigenval- 
ues by truncating the expansion for dct[7Y — E\ at 0{u^). 

Just as the dispersion relation for low lying states could 
have been derived by setting det[H — E] = to 0{v?) 
and solving the resultant quadratic equation for the low 
lying energy states: k = q-\-5 with small 5. Here we may 
do the same. 

The characteristic equation is easily read off: 

= det[H -E]= n^i[(w(fc,) - Ef\ - 

Au'^aj Y^[{v{kj,) + 8u){v{kj, + 2qi) + 8u)] 
h 

-Au^a\ Yy<^^h) + 8w)(w(% + 2^2) + 8u)] 
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x^^^^,{(v(k,)-Ef\ (130) 



where wc have shifted by a constant. Notice that 
decoupling trivially occurs - terms of the form [^^(A;) — 
Ey^ik!) — E\ where the modes k — k' = qi ± q2, cancel. 

For low lying states i.e. for the terms containing ki = 
qi + Si and ^2 = ^2 + S2: 



E = E„ 



+ a?[A||<52^i+Ax<5l;i] 
+4[A\\6l.2 + A±Sl;2] 



(131) 



with (5||,_L:m parallel and perpendicular to the minimizing 
manifold M at g^, trivially satisfies det[H — E] =0. 

This dispersion relation agrees, once again, with the 
result derived from the spin wave stiffness analysis 



^H=^J2(^{k)-v{q))\S{k)\' 
k 



(132) 



and by expansion of AH for different sorts of twists, the 
dispersion relations of the two spiral simply lumped to- 
gether. When d = 3, as in the 0{2) case {p = 1) this 



dispersion gives rise (in the Gaussian approximation) to 
diverging logarithmic fluctuations: 0(|ln e|). Apply- 
ing equipartition, the Gaussian spin fluctuations in the 
[2i— 1, 2i\ plane: 



A^2i-l,2i] = 0) > ASl^ pi-l,2i] = 0) 



= kHT 



(Pk 



(27r)3 a2[A||52.+A^5i^.]- 



(133) 



For all odd n and for all even n with p < n/2 there will 
be divergent fluctuations similar to those encountered for 
the 0(3) model. 

We can now update and summarize our conclusions: 
We have just proved that if frustrating interactions cause 
the ground states to be modulated then the associated 
ground state degeneracy (for n > 2) is much larger by 
comparison to the usual ferromagnetic ground states. For 
even n we have found that, generically, the a three dimen- 
sional system will not have long range order when M is 
two dimensional. When n is odd the system will never 
show long range order if M is (d-1) or (d-2) dimensional. 

If. in the continuum limit, the uniform (ferromagnetic) 
state is higher in energy than any other state then, by 
rotational symmetry, the manifold of minimizing modes 
in Fourier space is (d — 1) dimensional. 

In reality, small symmetry breaking terms (e.g. A 7^ 
in VQ{k)) will always be present- these will favor order- 
ing at a discrete set of {±9to}|^i- It n > 2\M\ then, 
irrespective of the even/odd parity of n, then will be a 



(AS'^ix^O)) > {n-2\M\) 



dPk 



knT 



{2t^Y v{k) - v{q) 



■ (134) 



This contribution is monotonically increasing in n; 
Within our scheme, Tc is flnite and may be estimated 
by the temperature at which the fluctuations are of or- 
der unity. By tweaking the symmetry breaking terms to 
smaller and smaller values, the fluctuation integral be- 
comes larger and larger. For instance, if take A 4; 1 in 
VQ{k) then the integral is very large and Tc extremely 
low (in can be made arbitrarily low). Thus as the system 
will be cooled from high temperatures, it might flrst un- 
dergo a Kostcrlitz-Thoulcss like transition at Tkt to an 
algebraically ordered state and develop true long range 
order at critical temperatures Tc < Tkt- 



XIII. LARGE N LIMIT 

So far we have seen that the even n systems are more 
"gapped" than its odd counterparts. There is never a 
paradox in the large n (or spherical model) limit. In this 
limit wherein a single normalization constraint is imposed 



(135) 
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the effective number of spin components n is of the order 
of the number of sites in the system N. The span of the 
system N (the number of Fourier modes allowed within 
the Brillouin zone) is always larger than the number of 
minimizing modes {qi}- 

In such a case we will be left with a divergence as in 



equation (134) due to the many unpaired spin compo- 
nents. 

In fact, within the spherical model, which is easily solv- 
able the fluctuation integral exactly marks the value of 
the inverse critical temperature 



1 



1 



.z. (STr)"^ vik)-v{q)' 



(136) 



Thus, Tc = if the latter integral diverges and our 
circle of ideas nicely closes on itself. 



XIV. 0(N > 2) WEISS MEAN FIELD THEORY OF 
ANY TRANSLATIONAL INVARIANT THEORY. 



Let us begin by examining the situation simple spiral 
states. In this case for 0{n > 2) when T < T^. 

{S{x)) = S Sa^ound-statei^^y (^37^ 

For the particular case 

c^qround— state / ->\ / -^ -^\ 

bl [x) — cos(q • x) 

aground- state,-.-. _ ■ 

ciQroujid— state / -t\ rt /i oo\ 

S!>3 [x) = (138) 

Now only the zbg modes have finite weight. Repeating 
the previous steps 

J2 V{x - 0, = (139) 



\{S{x^O))\ = \{S,{x = 0))\ 



(140) 



Define 



M[z] ^ ln[(2/z)("/2-i)/„/2_i(z)], (141) 

with [In/2-i{z)] a Bessel function. The mean-field equa- 
tion reads 

\{S,{x^O))\ = s = M[\Y,V{x-, y) {S, (y)) \ ] . 



The onset of the non-zero solutions is at 

\f3cvm=n. (142) 
li V{x — 0) — Q (no on-site interaction), then 



d-^k v{k) = 0, 



(143) 



implying that v{q) <Q and Tc > 0. Note that within the 
mean field approximation, Tc is a continuous function of 
the parameters. 

Here the ground state is symmetric with respect to all 
sites. The above is the exact value of Tc within Weiss 
mean field theory for the helical ground-states. 

For poly-spirals we will get p identical equations: both 
sides of the self consistency equations are multiplied by 
of where ai is the amplitude of the I — th spiral in the 
[(2/ — 1), 21] plane. As v{q,ri) — v{q), we will arrive at the 
same value of Tc as for the case of simple spirals. 



XV. EXTENSIONS TO ARBITRARY TWO SPIN 
INTERACTIONS: SPIN GLASSES, \U) SPACE 
(g)0{N) TOPOLOGY ETC. 

Any real kernel V{x, y) may be symmetrized \y{x, y + 
V{y,x)\/2 ^ V{x, y) to a hermitian form. 

Consequently, by a unitary transformation, it will be- 
come diagonal. The Fourier modes are the eigen-modes of 
V when it is translationally invariant. We may similarly 
envisage extensions to other, arbitrary, V{x, y) which will 
become diagonal in some other complete orthogonal basis 
\u): 



{ui\V\uj) = 5ij{ui\V\ui) 



(144) 



Many of the statements that we have made hitherto 
have a similar flavor in this more general case. 

For instance, the large n fluctuation integrals are of 
the same form 



d'^u 



1 



v{u) - Vn 



(145) 



with the wave-vector k traded in for u. 

Once again, one may examine the topology of the min- 
imizing manifold in u space. If the surface if (d — 1) di- 
mensional and v{u) is analytic in its environs then, for 
large n, Tc = 0. 

The topology of the ground state sector of 0{n) models 
will once again be governed by a direct product of the 
topology of the minimizing manifold in u space with the 
spherical manifold of the 0{n) group. In the general case 
in will be dramatically rich. 

We may similarly extend the Peierls bounds to some 
inflnite range interactions also in this case by contrasting 
the energy penalties in the now diagonalizing u basis with 
those that occur for short range systems in Fourier space 



^Edisordered = ^ I "^C^) 1^ ( ("I 



IS I "^Z ^mzn 



>^Y.\^^^)\>^^^ort{k)-V, 



short mini 



(146) 



18 



if they share the same lowest energy eigenstate of V. 



XVI. 0(N) SPIN DYNAMICS AND 
SIMULATIONS 



For our repeated general Hamiltonian 

H=lj2^{S,y)S{S).S{y} (147) 



the force on given spin 



dH 

' dS{z) 



For V^(.T, y) = V{x — y){= V(ij — x) or otherwise we may 
symmetrize \y{x,y + V{y,x)\/2 V{x,y), as we have 
done repeatedly, without changing H). For a more gen- 
eral two spin kernel V{x,fi) which is not translationally 
invariant the Fourier space index k should be replaced by 
the more general u. 



= _J2vix-y)S{y) 



d^Sjx) 

d^Sjk) 
dt^ 



= --v{k)S{k). 



(149) 



Alternatively, the last equation can be derived by starting 
with the Hamiltonian expressed directly in Fourier space 



df^ 



k 

dU{p,t) 

dH 



dt 



dS{p) 2iV 



[v{p) + v{-p)]Sip) 



= --v{p)S{p), (150) 



where in the last equation the momentum n(^) conjugate 
to S{x) is trivially 

dL 



d{dS{x)/dt) 

d{\ i:s[dS{x)/dtf \ V{x, y)S{x) ■ S{y)) 



d{dS{x)/dt) 



dSjx) 
dt 



, (151) 



and upon Fourier transforming n(p) = {dS{p)/dt). Let 
an arbitrary A satisfy 



A> — m.m{v{k)}. 

k 



(152) 



The equations of motion 
d'^S{k) 



^ = -[A + v{k)]S{k) 



(153) 



may be trivially integrated. [Adding the constant A 
merely shifts H ^ H + A/2.] 

^ ^ ^ ^ dSCk, t) 
Sun{k,t) = S {k, 0) cos uJkt-\ — |t=o X uj^^smujkt, 

Sunik, t + 6t) = S{k, t) cos ujkSt 
6S{k,t) _i . 



+- 



St 



This suggests the following simple algorithm: 

(i) At time t, start off with initial values {S{x,t)}. 

(ii) Fourier transform to find {S{k,t)}. 

(iii) Integrate to find the un-normalized {Sun{k,t + 

St)}. 

(iv) Fourier transform back to find the un-normalized 
real-space spins {Sun{x,t + 6t)}. 



(v) Normalize the spins: 

Six,t + 6t) = i^* + ^*^ 



\Sun{x,t + 6t)\ 



(154) 



(vi) Compute {6S{x,t+5t)} = {[S{x,t+5t)-S{x,t)]}. 

(vii) Fourier transform to find {6S{k, t + 6t)}. 

(viii) Go back to (ii). 



Thus far we have neglected thermal effects. To take 
these into account, one could integrate these equations 
with a thermal noise term augmented to the restoring 
force 

= -^v{k)S{k) + F^-^-(T,i). (155) 



Expressed in this format, the execution of this algo- 
rithm for continuous 0{n > 2) spins seems easier than 
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that for a discrete Ising system. Here the equations of 
motion may be integrated to produce arbitrarily small 
updates at all sites. 

This could, perhaps, be better than a brute force ap- 
proach whereby the torque equations in angular variables 
(the spins S{x) are automatically normalized) are inte- 
grated whereby the eqs of motion would explicitly read 



^V{x-y)sm[(j){x) - (/-(j/)]. 



(156) 



for an 0{2) system. 

For the three-component spin system: 



sm ti(x} 



sin 20{x 



d<l>{x) de{x) 



dt^ ' ' dt dt 

V{x — y) sin 9{x) x sin 9{y) sn\[(j){x) — 0(y)]; 



d'^e{x) _ 1 
~dtP~ ~ 2 



sin2.(f)(^)^ 

+ ^ V{x - 27){sin 6'(f ) 
y 

cos0{y) — cos9{x)sm9{y)cos[(f>{x) — ^(y)]}- (157) 



This may be computed by first integrating over fc|| em- 
ploying 



1 ~ cos\a(b ~ x)] TT , , 

^2 + ^2 ^^ rfx = -[l-e— cos(a6)] (162) 



to obtain 



A 



II -'0 



1 — cos(fc||a;|| + k± ■ x±) 

fc2 + (A^fclM||) 



dkw 



2k^\lAA Ll-exp(~J^A:ia;||)cos(A:'_L -fii)] (163) 



1 TT / 1 

"J. 



V ^1 

If (f) denotes the angle between k± and xj_ then 
l~A±^ 

1 — exp(— . / — — fcj^xii) cos(fc^ • xj_)]d(f> 



II 

\ A\^ 

27r — cxp(— , / — — /cj^xii ) / cos(fc^x± cos 0)c?(/). (164) 
>' A\\ > 



As 



Jo(x) ~ — cos(xcos0)d0. 



(165) 



XVII. APPENDIX 

Here we follow the beautiful treatment of Als- Nielsen 
et al. 0- 

Within the (hard spin) fully constrained XY model: 
G{x -y)^ {S{x) ■ sm = (cos[0(f ) - 9{y)])- (158) 

Here 9{x) ~ q ■ x + A9{x), i.e. A9 denotes the phase 
fluctuations about our spiral ground state and 

G(f _y) =cos(q- (f-y)))(e'(^«(^)-^«(^^)). (159) 

In our harmonic approximation {69{x)} are random 
Gaussian variables and only the first term in the cumu- 
lant expansion is non- vanishing. 
The correlator 

G{x) ^ cxp[-^[{\A9{x) - A9{0m 



exp 



keT 



d'^k 



1 — cosq ■ X 



{2ttY A\\5l^A^5\ 



(160) 



Now let us shift variables k ^ k — q = S, and for purposes 
of convergence explicitly introduce an upper bound on 

k±: < fc_L < A 



/(f) 



2Ai 



-7r27r 



A 1 — exp(— ^/^fc^a;||) Jo(fcj_a;^) 



-fc I dk \ 



We may now insert the series expansion of Jo (a;) and 
integrate term by term. 



(-)"z2 



22"(n!)2 



(166) 



Comparing the result to the series form for the exponen- 
tial 



E^{z) = -^-\nz-Y,{' 



n{n\) 



(7 is Euler's constant) we find that 

Ad^ n 
G(f)^ — exph2,77-^i?i( ^ 



(167) 



cos[(7X||] (168) 



where ri = ^^\/4^ which in our case is -If^Q^^^ and 

d — ^ where A is the aforementioned ultra violet mo- 
mentum cutoff. 



I{x±,xi\) 



1 — cos{q ■ x) 



(161) 
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